ICM 2002 • Vol. Ill • 1-3 



The Dynamical Systems Approach to the 
Equations of a Linearly Viscous 
Compressible Barotropic Fluid 



E. Feireisl* 



Abstract 

We develop a dynamical systems theory for the compressible Navier-Stokes 
equations based on global in time weak solutions. The following questions will 
be addressed: 

• Global existence and critical values of the adiabatic constant; 

• dissipativity in the sense of Levinson - bounded absorbing sets; 

• asymptotic compactness; 

• the long-time behaviour and attractors. 
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1. Introduction 

The long-time behaviour of solutions to the evolutionary equations arising in 
the mathematical fluid mechanics has been the subject of many theoretical studies. 
This type of problems is apparently related to the phenomena of turbulence, and 
there is still a significant gap between many formal "scenarios" and mathematically 
rigorous results. 

The dynamical systems in question are usually related to a system of partial 
differential equations and, consequently, they are defined in an infinite dimensional 
phase space. On the other hand, the presence of dissipative terms in the equations 
due to viscosity results in the existence of global attractors-compact invariant sets 
attracting uniformly in time all trajectories emanating from a given bounded set. 
Such a theory is well developed for the incompressible linearly viscous fluids, and 
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the reader may consult the monographs of BABIN and VISHIK [T], TEMAM [211 
or CONSTANTIN et al. [2] for the recent state of art. 

On the other hand, much less seems to be known for the compressible flu- 
ids. While there is an existence theory of the weak solutions for the incompressible 
Navier-Stokes equations due to LERAY ^Hli its "compressible" counterpart ap- 
peared only recently in the work of LIONS jJUl- Even in the incompressible case, 
there is a qualitative difference between the two-dimensional case solved by LA- 
DYZHENSKAYA jEl , and the three-dimensional case representing one of the most 
challenging unsolved problems of modern mathematics. It is worth-noting that a 
similar gap divides the one and more-dimensional problems for the compressible 
fluids. 

The time evolution of the fluid density g = g{t, x) and the velocity u — u{t, x) 
is governed by the Navier-Stokes system of equations: 

5te + div(ei?) = 0, (1.1) 

dt{Qu) + div(£»u ® u) + Vp = div S" + gf, (1.2) 

where p is the pressure, S the viscous stress tensor, and / a given external force. 

In what follows, we consider linearly viscous (Newtonian) fluids where the 
viscous stress is related to the velocity by the constitutive law 

S = II {Vu + Vu^) + A div u J, (1.3) 

where the viscosity coefficients satisfy 

/i > 0, \ + ^i>0. (1.4) 

Generally speaking, the pressure p depends on the density and the internal en- 
ergy (temperature) of the fluid. If it is the case, the system l|l.l|) . (|1.2|) is not closed 
and should be complemented by the energy equation. Unfortunately, however, the 
available global existence results for this full system allow for only for small initial 
data (cf. MATSUMURA and NISHIDA |1]]). 

There are physically relevant situations when one can assume the flow is 
barotropic, i.e., the pressure depends solely on the density. This is the case when 
either the temperature (the isothermal case) or the entropy (the isentropic case) 
are supposed to be constant. The typical constitutive relation between the pressure 
and the density then reads 

p = p{g) = ag'^, a > 0, (1.5) 

where 7 = 1 in the isothermal case, and 7 > 1 represents the adiabatic constant 
in the isentropic regime. More general and even non-monotone pressure-density 
constitutive laws arise in nuclear plasma physics (see @]). In the barotropic regime, 
the equations Hl.l|l . (|1.2|) form a closed system and complemented by suitable initial 
and boundary conditions represent a (at least formally) well-posed problem. 
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If the problem is posed on a spatial domain fl C R , one usually assumes that 
the fluid adheres completely to the boundary which is mathematically expressed by 
the no-slip boundary conditions for the velocity: 

u\an = 0. (1.6) 

Note that for viscous fluids such a condition is in a very good agreement with 
physical experiment. 

In accordance with the deterministic principle, the state of the system at any 
time t > to should be given by the initial conditions 

9{to) = Qi, {9u){to) = qi- (1.7) 

The function gj is non-negative and the momentum qj satisfies the compatibility 
condition 

qi = a. a. on the set {gi — 0}. (1.8) 

The reason why we impose the initial conditions for the momentum gu rather than 
for the velocity w is that the former quantity is weakly continuous with respect 
to time while the instantaneous values of the velocity are determined only almost 
anywhere with respect to time. Clearly such a problem does not arise provided the 
initial density gj is strictly positive. However, it is an interesting open problem 
whether or not this property is preserved at any positive time for any distributional 
solution of the problem satisfying the natural energy estimates (cf. HOFF and 
SMOLLER TO). 



2. Finite energy weak solutions and well-posedness 

In order to study the long-time behaviour, one should first make sure that the 
class of objects one deals with is not void. More precisely, one should be able to 
prove the existence of global-in-time solution for any initial data gj, qj satisfying 
some physically relevant hypothesis. 

Multiplying the equations of motion by u and integrating by parts one deduces 
the energy inequality 

^E[g,u\ + I ^|Vup + (A + ^)|div up dx < / gf-uAx, (2.1) 

where the total energy E is given by the formula 

E[g,u]^ / g\u\^ + P{g) dx 
Jn 

with 

P\g)g-P{g) ^p{g). 

Note that in the most common isentropic case, the function P can be taken in 
the form 
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The energy inequality is the main (and almost the only one) source of a priori 
estimates. Accordingly, "reasonble" solutions of the problem (II. 111 . (|1.2I) defined on 
a bounded time interval I C R should belong to the class 

g>0, ge L^iI;L'<{n)), u £ L''{I;Wo\n, R^)). (2.2) 

The energy inequality H2.1|l represents an additional constraint imposed on 
any solution g, u of the problem. Similarly as in the theory of the variational 
(weak) solutions developed for the incompressible case by Leray, it is not clear if it 
is satisfied for any weak solution of the problem. 

Following DiPERNA and LIONS |H] we shall say that m is a renormalized 
solution of the continuity equation if the identity 

dth{g) + div(6(g)w) + {b'{g)g - b{g)j div u = (2.3) 

holds in the sense of distributions for any function b G C^{R) such that 

6'(^?) = for all > const(6). (2.4) 

Similarly as for the energy inequality, it is not known if any weak solution u of 
satisfies H2.3II . 

Definition 2.1 We shall say that g, u is a finite energy weak solution of the 
problem - I J. 6]) on a set I x fl if the following conditions are satisfied: 

• The functions g, u belong to the function spaces determined in 1^2. 

• the energy inequality 1^2.1]) is satified in "D'^I); 

• the continuity equation as well as its renormalized version \2.'^) hold in 
■D'{I X R^) provided g, u were extended to be zero outside Vt; 

• the momentum equation 1^1. 2^) is satisfied in T>'{I x f2). 

The most general available existence result reads as follows: 
Theorem 2.1 Let $7 C R^ , N = 2,3 be a bounded Lipschitz domain. Let 
I = (0,T), and let the initial data gi, qi satisfy ^.8\) together with 

gi£L\Vt), ^eLi(r!). 
Qi 

Let f be a bounded measurable function oftCzI,xCzQ. Finally, let the pressure 
p G C[0, oo) n C^{0, oo) be given by a constitutive law 

p — pig), —g^^^ — b < p'{g) < ag^^^ + b, for all g > Q, 
a 

where a > 0, 6 > 0, and 

N 

Then the problem admits a finite energy weak solution g, u on I xQ. 

satisfying the initial conditions 

In his pioneering work, LIONS ^H] proved Theorem 2.1 for regular, p mono- 
tone, and 7 > I if TV = 2, and 7 > | for A*" = 3. The hypotheses concerning 7 were 
relaxed in the case of a general bounded domain treated in JI], and the 

hypothesis of monotonicity of the pressure removed in 0]. 
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3. Ultimate boundedness 

The first issue to be discussed when describing the long-time asymptotics of a 
given dynamical system is ultimate boundedness or dissipativity. This means there 
exists an absorbing set bounded in a suitable topology. Here "suitable topology" 
is of course that one induced by the total energy E. One of possible results in this 
direction is contained in the following theorem. 

Theorem 3.1 Let Q C i?^, = 2,3 be a bounded Lipschitz domain. Let f 
be a bounded measurable function such that 

ess sup \f{t, x)\ < F. 

Let the pressure p be given by the isentropic constitutive law 

5 

p = ag^ with ^ > I if N = 2, j > - for N = 3. 

o 

Finally, set 

/ g dx — m > 0. 
Jn 

Then there exists a constant E^o depending solely on m and F having the 
following property: 

Given Ej , there exists a time T = T{Ei) such that 

E[Q,u\{t) < Eoc for a.a.t > T 

provided 

esslimsupi?[p, 'it](t) < Ej 

t-*o+ 

and Q, u is a finite energy weak solution of the problem - \1.6]) on (0, oo) x fl. 

The proof of Theorem 3.1 can be found in and 7 . The reader will have 
noticed the "critical" exponent 7 > | which is larger than in the existence Theorem 
2.1 and, as a matter of fact, does not include any physically relevant case. Indeed the 
value I happens to be the largest adiabatic constant corresponding to a monoatomic 
gas. 

Under the hypotheses of Theorem 3.1, the dissipative mechanism induced by 
viscocity is strong enough to prevent any "resonance" phenomena, i.e., the existence 
of unbounded solutions driven by a bounded external force. Another interesting 
feature is the existence of periodic solutions provided / is periodic in time. 

Theorem 3.2 In addition to the hypotheses of Theorem 3.1, assume that f is 
time periodic, i.e., 

f{t + uj,x) = f(t, x) for a.a. t G R, x G 

with a certain period lu > 0. 

Then there exists at least one finite energy weak solution of the problem fl.l\ - 
\1.6]) on R X fl which is lu— periodic in time, i.e., 

g{t + uj) = g{t), {gu){t + uj) = {Qu){t) for all t G R. 



See ^U] for the proof. 
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4. Asymptotic compactness 

The property of asymptotic compactness of a given dynamical system plays 
a crucial role in the proof of existence of a global (compact) attractor. Here, the 
main problem is the density component which is bounded only in L'^{Q), and, con- 
sequently, compact only with respect to the weak topology on this space. Moreover, 
given the hyperbolic character of the continuity equation, one cannot hope any pos- 
sible oscillations of the density to be killed at a finite time. However, the amplitude 
of possible oscillations is decreasing in time uniformly on trajectories emanating 
from a given bounded set. 

Theorem 4.1 Let il C , N = 2,3 be a bounded Lipschitz domain. Let the 
pressure p be given by the isentropic constitutive relation 

5 

p = ag'', a > 0, 7 > 1 i/iV = 2, 7 > - /or iV = 3. 

Let fn be a sequence of functions such that 

ess sup |/n(i,2:)| < F independently of n ~ 1,2, ... . 
teR. xeii 

Finally, let Qn, Un « sequence of finite energy weak solutions to the problem 
- \1.6]} on (0,00) X fl such that 

Qn da; = m, 

ess lim sup i?[f3n, u„](<) < Ei 

independently of n = 1,2, ... . 

Then any sequence of times t„ — > 00 contains a subsequence such that 

Qn{tn + t) ^ g{t) strongly in L^{Q) and weakly in L'^ (Q), 

{QnUn){tn + t) {gu){t) wcakly in L^{Q,R^) 



for any t R. 

Moreover, 



\(>n(tn +t,x) — g(t,x)\'^ dxdt 



J Jn 



\{gnUn)itn +t,x) ~ {gu){t,x)\ dxdt 

ij Jn 

for any bounded interval J <Z R. The limit functions g, u represent a globally defined 
(for t 1^ R) finite energy weak solution of the problem \1.1\ - TT1^) driven by a force 



f — lim /„(t„ + ■) in the weak star topology of the space L°°{R X n,R'''), 

n — >oo 

and such that 



eas snp E[g,u](t) < 00. 
tefl 

The proof of Theorem 4.1 is given in ^] (see also 
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5. The long-time behaviour, attractors 

The results presented in Sect ions lTl ITI allow us to develop a theory of attractors 
analogous to that one for the incompressible flows (see e.g. TEMAM [21] )• Assume, 
for the sake of simplicity, that the driving force / is independent of t. We introduce 

"4 = {[gi, qi] I gi — £'(0), qi ~ {gu){0) where g, u is a finite energy weak solution 



of the problem mU-Ol) on R x il with E[g,u] e L°°{R)}. (5.1) 

In other words, the set A is formed by all globally defined (for t E R) trajec- 
tories whose energy is uniformly bounded. 

The next result shows that ^ is a global attractor in the sense of FOIAS and 
TEMAM 

Theorem 5.1 Let O C R^ , N — 2,3 be a bounded Lipschitz domain. Let 
f = f{x) be a bounded measurable function independent o J time, and let the pressure 
p be given by the isentropic constitutive law 

p = ag'^, a > 0, -f > 1 for N = 2, 7 > | if N = 3. 

Then the set A defined by 1(5.1}) is compact in the space 

and 

sup inf (\\git)-gi\\L''{n)+ iigu){t) - qj) ■ (/) dx ) 

for i — > oo 

for any 1 < a < 7 and any (f> G L^-i (fi, R ), where the symbol B{Ej) stands for 
the set of all finite energy weak solutions of the problem 1) j. Jl - li.6|) on (0, 00) x Q 
such that 

ess \im sup E[g,u]{t) < Ej. 
The proof can be found in [S]. 

To conclude, we give another result concerning the long-time behaviour of 
solutions on condition that the driving force is a gradient of a scalar potential. 

Theorem 5.2 Let il C R^ , N — 2,3 be a bounded Lipschitz domain. Let the 
pressure p be given by the isentropic constitutive relation 

p = ag^, a > 0, 7 > y- 

Let the driving force f be of the form 

f = VF, 
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where F = F(x) be a scalar potential which is globally Lipschitz on f7 and such that 
the upper level sets 

[F>k]^{xC,n\ F{x) > k} 

are connected for any k £ R. 

Then any finite energy weak solution g, u of the problem 1^1. 1\ - 1,?. 6)) satisfies 

g{t) Qs in L'^{n), {gu){t) in L^[Q) as t ^ oo, 

where Qs is a solution of the static problem 

ay q1 — Qsf on VL. 
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